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Abstract 

The scattering problem under the influence of the Aharonov-Bohm (AB) poten- 
tial is reconsidered. By solving the Lippmann-Schwinger (LS) equation we obtain 
the wave function of the scattering state in this system. In spite of working with a 
plane wave as an incident wave we obtain the same wave function as was given by 
Aharonov and Bohm. Another method to solve the scattering problem is given by 
making use of a modified version of Gordon's idea which was invented to consider 
the scattering by the Coulomb potential. These two methods give the same result, 
which guarantees the validity of taking an incident plane wave as usual to make an 
analysis of this scattering problem. The scattering problem by a solenoid of finite 
radius is also discussed, and we find that the vector potential of the solenoid affects 
the charged particles even when the magnitude of the flux is an odd integer as well 
as noninteger. It is shown that the unitarity of the S matrix holds provided that a 
plane wave is taken to be an incident one. 
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1 Introduction 



Since Aharonov and Bohm (AB) have discussed a scattering problem of a charged par- 
ticle by a solenoid in order to clarify the significance of the vector potential in the quantum 
theory many people have considered the same problem from various viewpoints @- 
|19H . As is well known, there are two approaches to deal with a scattering problem in the 



quantum theory. The first approach is to find a stationary state describing the scattering 
process by solving a time-independent Schrodinger equation. The second one is to study 
the time development of a wave packet with respect to a time-dependent Schrodinger 
equation. Most people as well as Aharonov and Bohm have analyzed the scattering by 
means of the first approach, and some people have discussed the same problem with the 



second approach pT| , [[i9|j . As we see in the following, however, in spite of these efforts it 
seems not to be clear what is an incident wave in this scattering process. In this paper we 
will try to answer this question with the first method because there seems to be lacking 
a common interpretation of the stationary wave function in the literature. 

A system of charged particles interacting with the solenoid is described by a Hamil- 
tonian 

where the electromagnetic vector potential A is given by 

$ 

A(x) = (— sin (p, cosy?) , (x, y) = (r cos cp, r sin cp). 

2nr 

In order to study the scattering of the charged particles we solve the time-independent 
Schrodinger equation 

HMr,<p) = Eil> E (r,<p) (1-2) 
to find an eigenfunction which describes the scattering process of charged particles. 

Since the Hamiltonian commutes with the angular momentum, it can be easily shown 
that a most general solution of (pL.2|) is given by 



Mr,<p)= E c n e in *J ]n+al (kr), E=^-, (1.3) 



-co 



where J u (x) denotes the Bessel function of u-th order and we have put a = —e^/2-K%c. 
In ( p.. 3D , c„'s are arbitrary constants to be determined by phisycal requirements. Our 
main interest here and in the following is to ask what should be the correct choice for the 
coefficients c n 's to describe the scattering process. 

It has been asserted by Aharonov and Bohm and many other authors that the incident 
wave for this scattering problem, when the incident beam comes into from the positive 
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x axis, should be a modulated plane wave e lkrcost f ioll p to make the probability current 
of the incident wave constant. To fulfill this requirement the coefficients have been taken 
to be c n = {-i)\ n+a \. In the case of nonintegral a, however, the incident wave becomes 
a multi-valued function. We may say that from the physical point of view it seems 
quite unsatisfactory to take such a multi-valued wave function as an incident one. We 
should also note that, for sufficiently large r, the term due to the vector potential do not 
contribute to the dominant part of the current, even if we take a plane wave as an incident 
wave function. Thus, it is not conclusive to argue that the incident modulated wave gives 
the condition to determine those constants as c n = (— i)' n+a L 

It will be, therefore, instructive to reconsider the scattering problem from other view- 
points. In this paper we try to find the wave function to describe the scattering process 
using two different ways. The first one is to solve the Lippmann-Schwinger (LS) equation, 
which will be a standard method to consider the scattering problem of a quantum system, 
taking a plane wave as an incident state instead of the modulated one. If we adopt the 
Born expansion to solve the LS equation, we will soon meet a difficulty because the per- 
turbative method does not work for the present problem as is shown in fll5| , Then we 
find an exact solution of the LS equation with the aid of the Feynman kernel. The second 
method is an application of Gordon's idea which is proposed to discuss the scattering 
problem by the Coulomb potential |13|| . The idea may have been introduced to avoid 
the difficulty caused by the long-ranged nature of the Coulomb potential in formulating 
a scattering theory. This method will also be useful to examine scattering problems by 
other potentials with long-range effect. It will be shown that these two methods give the 
same wave function to describe the scattering state given by AB [|I|]. 

If a = integer, a further discussion is needed since the solution of the Schrodinger 
equation for infinitely thin solenoid does neither vanish nor be defined at the origin in 
that case. To ensure the impenetrability of the solenoid even for integral a, we assume a 
solenoid to have a finite radius and consider this issue by generalizing the second method. 

The plan of this paper is as follows. In section ^| we give the Feynman kernel with 
effects of the solenoid potential using the path integral method. Section |3| is devoted to 
solving the LS equation exactly with the aid of the Feynman kernel, and the explicit form 
of the wave function will be obtained. In section |] Gordon's method will be argued and 
its generalization to a system of a solenoid with finite radius will be done in section [5| 
Conclusions and discussions in comparison with other's results are made in section 0. 
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2 The Feynman kernel with effects of the solenoid potential 



Since the complete set of the eigenfunctions 

1 



'2ix 



e m ^J ln+al (kr), n = 0,±l,±2,... 



for the Hamiltonian with effects of the solenoid are known, we can immediately find an 
expression of the Feynman kernel 

K(x F , xj; T) = (x F \e- lflT/h \ Xl ) (2.1) 

as its spectral representation |§, |ll 



roo 1 +°° 

K(x F ,x r ,T)= k dke- ink2T/2 »— Yl e in ^-^J ln+a \(kr F )J ]n+al (krj). (2.2) 

J 27r n=-oo 

To verify this expression, it suffices to notice the following facts, (i) it obeys the time- 
dependent Schrodinger equation, (ii) it is apparently single-valued with respect to both 
x F and Xj, (iii) it has the correct limit 

\im K(x F , Xj;t) = 5 2 (x F — Xj), (2.3) 



1 r°° 

—S( a -b)= kdkJJak)JJbk) \Re(u) > -1, a,b > , (2.4) 
ab Jo 



which follows from 
1 

\fab 
and 

+oo 

e in0 = 2n5(9) [-tt < 9 < it) . (2.5) 

71= — CO 

If we carry out the integration with respect to k in ( p.2[) , we obtain 



n=— oo 



By use of (|2.2|) or (|2.6|) we can proceed to solve the LS equation. In this section, 
however, we would like to give another derivation of (|2.2|) because there seems to be some 
confusion in the path integral construction of the Feynman kernel in the literature J7|- 
||10|| . Among them the most typical one would be the interpretation of its expression as 
the sum over winding number such as 



K(x F ,xj;T)= Y e- ia ^-^- 2m ^K m (x F ,x i; T), (2.7) 



m=— oo 
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where K m (xF, Xf, T) comes from paths going around the solenoid m times in anticlockwise 
way, and the factor e 2mma is a one dimensional representation of the fundamental group 
of the configuration space. In the following, however, we show that the sum over winding 
number is not essential, and that the sum in ( J2.7| ) is to be interpreted as a result of the 
reformulation with aid of the Poisson sum formula for the expression obtained using the 
usual path integral method. To achieve this we will make use of the completeness relations 
of both eigenvectors \x) of x with eigenvalue x, and eigenvectors \p) of p with eigenvalue 
p, in formulating the path integral. As is easily recognized, the single-valuedness of (|2.7| ) 
is the consequence of these relations [12[ . 



Now we give the path integral derivation for the Feynman kernel of this system. Defin- 
ing the exponential operator by 



-iHT/h 



lim (1--H 



N 



T/N, 



(2. 



and using the completeness of the states \x), we obtain 



r ( is - \ 

K(x F , Xl ;T) = JimJ [] E|>0') I (l - jHj \x(j - 1)). (2.9) 



Then, with the aid of the completeness of the states \p), we can express the infinitesimal 
version of the Feynman kernel as 



(x 



lim 



d 2 p 



exp < -p(x 



x 



(27Th) 2 

where A(x, x') = {A(x) + A(x')}/2. 

After shifting the integration variable p by 



P 



is 1 



ip- 



, (2.10) 



e— r, 



p i— >■ p H — A(je, cc' 



we can rewrite (|2.10|) as 

is 



(x 



h 



H)\x') = exp < —A(x,x')(x - x' 



he 



x lim 

<5^0 



d 2 p 
(2nh) 



% 

h' 



2 exp <; -p{x - x 



i KM (2 - n; 



Carrying out the Gaussian integration with respect to p and noting that 



A(x, x')(x - x') 



$ fr' r\ . , 



(2.12) 
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we obtain 



is 



(x\{l--H)\x') 



fie iS 
hm — exp 

5^0 2mne 



~2he~ 



■{r 2 + r' 2 — 2rr' cos(f — </?')} 



la I r r \ . . . 



(2.13) 



where we have made a change of variables from Cartesian coordinate to the polar one and 
5 in (|2.13|) has been renamed by 1 — ihpS/e >— > e~ lS . 

To find the kernel for a finite time interval T we need to perform N — 1 integrations 
with respect to £c's in (|2.9|) . For this purpose, the form of the exponent in (|2.13|) is 
extremely inconvenient. To overcome the difficulty it is useful to rewrite 



he 



V rr c i8 

he 



a 


( r' 


' 2 


. r 


- 4 


r 


r 





1/2 



cos(v? -if' - 9 £ ) 



^Vyl + tan 2 9 e cos(<^ -if'- 9 e ), 



where 



tan 6L 



hea 



-id 



(2.14) 
(2.15) 



Upon integration with respect to x or x', the Gaussian part in the integrand will dominate 
for sufficiently small e. We may, therefore, regard components of X X clS 0(e 1/2 ). Then 
it follows that 

(2.16) 



U- + L )=l + 0(e^). 



2 \ r r 1 

(Note, however, that the same argument does not hold true for if — if'.) Recalling that we 
may discard terms of 0(e p ) for p > 1, in the exponent of a path integral, we can replace 
the definition of tan 9 £ by 

tan £ = ^ e - i5 f 1 + O^e 1 ' 2 )) = 9 £ {l + 0(e^ 2 )} . 



Thus we obtain 



exp 
exp 



•A 477 " i<5 / /\ tat r r \ 
~ l ~he~ e C ° S ^ _ ^ ' ~ ~2 [ 7 V' I ~ V ' 



-i^-e i5 Jl + 9 2 e cos((p - f' 



he 



{l + 0(e 3 / 2 )} 



+oo 



\n\ 



ihe 



e i5 sfl + 9 2 e *»(v-v'-*«) {l + 0{e 3/2 )} . 



(2.17) 
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When e becomes small, the argument of the modified Bessel function grows to allow us 
to apply its asymptotic form. Then keeping terms up to 0(e) in the exponent, we obtain 



>| 



ihe 



e~ lS exp 



2irnrr' 

' \trr j S 



firr iS 



ihe 



e — 



-^{(n + af-l/A) 



ihe 2/irr 
e iS ) {l + 0(e 3 / 2 )}. 



{l + 0(e i ' 2 )} 



Substituting it into (|2.17| ), we arrive at 



(2.18) 



cxp 



— i— — e cos(<^ — ip ) 



la, I r 



he 



+ - sm(tp - <p') 
2 \ r r 



>^'){l + 0(£ 3 / 2 )} . 



(2.19) 



Therefore the infinitesimal kernel ( |2.13| ) is now rewritten as 



K(x, x'\ e) 



le 



(x\[l--H)\x') 



fie 



lim 

2-Kihe 



53 ex p 



/le 1 



2ife 



(r 2 + r' 2 ) + m(<£ - y?') 



xJ, 



n+a| 



ihe 



e i5 ){l + 0(e 3 / 2 )} 



(2.20) 



Here a comment is in need; in obtaining the result of ( |2.18j ) we have discarded the possi- 
bility to use the modified Bessel function of negative order since it breaks the regularity 
of the kernel at the origin. 



(2.21^ 



It is now straightforward to see that the multiplication rule holds: 

J d 2 xK(x2, x; e)K(x, X\;e) = K(x 2 , X\\ 2e), 

since the integration with respect to the angle variable is trivial and we may make use of 
a formula 

/ rdre- ar J r >r) J ^qr) = —e^ +c > )/ 4n / /t (p 9 /2o), (2.22) 
jo 2a 

which holds for |arg(a)| < 7r/2, Re(/i) > —1, p, q > 0. Repeated use of the rule ( |2.21j) 
(and putting all S's to after integration) will lead us to 



K(x,x';T) 



2nihT 



53 exp 



f in 
\2hT 



(r' + r' 2 )\l ln+al 



ihT 



(2.23) 



Thus (|2.6| ) is again obtained by the usual formulation of path integral. Here it should 
be noticed that the sum over winding numbers in formulating the path integral is not 
essential. 
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3 The wave function of a scattering state as a solution of Lippmann- 
Schwinger equation 



In this section we obtain the wave function for the scattering state of charged par- 
ticles scattered by the solenoid. It is known that, for the present problem, the Born 
approximation fails to give a reliable answer because we cannot avoid a divergent integral 
Jq dr' J^kr') / r' even in its first order []15 , [16] . The iterative method to solve the LS 
equation will also be unsatisfactory by the same reason. Therefore we need to solve it in 
an exact way with the aid of the Feynman kernel given in (|2.23 ). 

The LS equation for the system reads 

if>E(r,ip) = u E (r,(p)+ip s (r,ip), (3.1) 

r°° f +7T / " \ — i h 2 

ip S (r,<p) = / r'dr' d<p' (x\ [E - H + is) \x') —^a (-2id v , + a) u E (r',(p') , 

JO J— 7T ZfiT 

where we have taken an incident plane wave Ue{t,<p) = e lhrcos8 (8 = ip — (fio) as an 
eigenstate of the free Hamiltonian and </2 indicates the direction of the incident beam. 
From ( J2.23[ ), we can easily obtain the Green's function in the above by Laplace transform 

(X\ (E-H+ teY 1 \x') = r ^L e We)T/K , l e -iHT/H ,/v ,3^ 

v ' Jo in 



By putting E = fr 2 fc 2 /2/i, it turns out to be 

JL 

2ih 



(x\(E-H + ie)~ \x') = E ^-^{8(r-r')H^ +al (kr)J ln+a] (kr') 



n=—oo 



I™ P dp \ {ap) f {bp) = ^Hl l \ak)J v {bk) \Rb{u) >-l,a>b> 0], (3.4) 

JO J) K 26 A 



+8(r'-r)J ln+al (kr)H^ +al (kr')}, (3.3) 

where use has been made of a formula 
J v {ap)J u {bp) _ iri 

and 8(x) is the step function. 

Substituting ( |3.3|) and partial wave expansion of the plane wave ue (r, <p) into the 
integrand of ipc (r, <p), we obtain 

^S V) = E { A n (0 KU ( kr ) + B - ( r ) J l-M (kr)}e md+ ^' 2 , (3.5) 

n=— oo 

where 



rr dr' 

A n {r) = —a(2n + a) —J\ n+a \{kr') J\ n \{kr') 



B n (r) -- 
nK ' 2i 



71 

2i y ' Jo r 

7T , r°° dr' 



/■OO H r ' , . 

a(2n + a) —H^kr 1 ) J ]n] (kr'). (3.6) 

J r T 
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Making use of a formula of indefinite integral for cylindrical functions(represented by Z^ 
and Z,j for the sake of convenience) 



/ ' —Zn{ax) Z u (ax) = aX {z^+i (ax) Z v (ax) - Z„ (ax) Z v+1 (ax)} 

J x u — v l 1 J 

(jm^v), (3.7) 



Z^ (ax) Z v (ax) 



we obtain 



fi + U 



ITT ( "i 

A n (r) = —kr {J\ n+a \ +1 (kr) J\ n \(kr) - J\ n+a \(kr) J\ n \ +1 (kr)j 



2 

in a (2n + a) 
2 \n + oil + In 



J\n+a\ (kr) J\ n \ (At) (3i 



and 

+TlT ! Sl<l| M J H M + e-^^l-"^ 2 . (3.9) 
2 |n + a | + |n| 11 

By a simple calculation with the aid of Lommel's formula 

J v+1 (x)H£\x) - J v (x)Hi%(x) = — , 

TlX 



we have 



[A n (r) H^ +al (kr) + B n (r) J\ n+a \ (kr)} e **+*M*/2 

- J w (Ax) e me+l ^ 12 + J\ n+a \ (At) e *"*+*^-*l«+°l*/2 . (3.10) 



Then ( |3.5| ) can be rewritten as 

+00 

^ s ( r , tp) = _ e ^co S + J2 J| n+a , (At) e W»^l™+°^/2. (3.11) 

n=— 00 

We thus find that the total wave function for the scattering state is given by 

+00 

iJ> E (r 1 <p) = £ J\ n+a \(kr)e me+in ^ n+a ^l 2 . (3.12) 

n=— 00 

It is very interesting to recognize that by putting ip = n the solution of LS equation 
coincides with the wave function obtained by AB Jl|. But we have to remember that in 
solving the LS equation we take the plane wave as an incident wave and that the resulting 
scattered wave is given by (|3.11|) . In spite of the fact that the total wave function is same 



as that of AB, both the incident wave and the scattered wave in this section are different 
from those of AB as a consequence. 
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Next we proceed to find the differential cross section by use of the scattered wave (|3.11|) . 
In view of Q and we notice that A n (r) = 0({kr) ) while B n (r) = 0((kr)- 1 ) for 

large kr. Then we easily obtain the asymptotic form of scattered wave (r, (p) from 
( p.5| ) as 

+00 

ifeMrZo S ^(oo)^ Q| (fcr)e^'^/ 2 , (3.13) 

n=— 00 

where A, (00) is found from (13.8!) to be 



A n (oo) = -isin{(|n + a\ - |n|)7r/2}. (3.14) 



Using the asymptotic form of the Hunkel functions and (|3.14j ) for A n (oo) in (|3.13|) , we 
are lead to 



1 +00 
(r, ip) ~ / 4 £ - l)e m9 , (3.15) 

V LlXkV n=-oo 

where the phase-shift in n-th partial wave is given by 

SM=l~^ (; + [«^°) . (3.16) 
V ; I +7ra/2 [n+[a] < 0) V ' 

Here and in the following we denote the integral part of a by [a] and its nonintegral part 
by {a} to write a = [a] + {a}. 

We here introduce a regularization parameter e for the sum in (B7T5) so that it is 



defined as an Abel sum because the phase-shift does not decrease at all when \n\ becomes 
large and define f(9) as 



g— i\a\t) ( +00 +00 



f(9) = lim-== I ^(e" mQ - l) e m(> - nt + ^(e*™ - i) e - me " ne . (3.17) 



vft=0 n=l 



Then performing the sum of geometric series and making use of a symbolic relation 

lim = P =F inSfx — a) 

e-*o x — a±ie x — a 

with denoting the principal value by P, we finally obtain in terms of the scattering am- 
plitude f(0) 

v r 

f(6) = ^{(cosna- 1)5(9) + i^e^P-^). (3.18) 

Although the total wave function has happened to have exactly same form as the result 
of AB as mentioned in the above, the scattering amplitude ( |3.18| ) disagrees to that of AB 



9 



0] or by the 5 function term. Nevertheless, this disagreement can be discarded when 
we have an interest in the differential cross section only for non-forward direction(# 7^ 0) 
since we cannot well separate the scattered and un-scattered particles in the forward 
direction experimentally. As far as in the non-forward direction, the differential cross 
section is thus given by 

da (a) = 1 Sm 9 ™j e . (3.19) 
v ; 2vrA;sin 2 (#/2) v ' 

In this sense the scattering amplitude of AB describes the physics appropriately. However, 
if we take into account the unitarity of the S'-matrix, the S function for the forward 
direction cannot be neglected as is pointed out by Ruijsenaars ||14|| . A more detailed 
description of the property of the ^-matrix for this system is given in appendix [FJ. 



4 Another derivation of the scattering state 



We here consider the other approach to the problem by using a modified version of 
Gordon's idea which has been proposed in the analysis for the scattering of a charged 
particle by the Coulomb potential ||13| . The essence of the method is to prepare the 
asymptotic region described by the free Hamiltonian far distant from the solenoid in 
order to overcome some difficulties caused by the long range effects of the solenoid field. 
To this aim we introduce a modified vector potential 




1 

_r_ 

R 







re 9 (0 < r < r ) 

e v {r < r < R) 
(R < r) 



(4.1) 



where r is the radius of the shielded solenoid. It should be noticed that in the region 
R < r the vector potential does not affect charged particles. To go back to the original 
AB problem we just put R — > 00 after solving the Schrodinger equation for this system. 
In this section, we first deal with the scattering by an infinitely thin solenoid(ro = 0), and 
then we generalize the analysis to the case of a finite size(r > 0) solenoid in the next 
section. 

In the asymptotic region(r > R) where the vector potential is absent, the solution of 
the Schrodinger equation is given by eigenstates of the free Hamiltonian and the wave 
function to describe the scattering state ipn(r, if) will be given by 



^11 ( r ,¥) = e 



ikr cos ( 



+ J2 OnHP(kr)e 



inO 



V 5 - <Po)i 



(4.2) 
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where a n 's are constant coefficients to be determined in the following. In the scattering 
region (0 < r < R) the wave function ^j(r, ip) is subjected to 

h 2 k 2 



Hi 



2fi 
2\i 



Assuming the partial wave expansion for ipi(r, ip) 



+00 



n + a I 1 — 



r 



fc 2 



r 







and making a change of variable r 



z = a{r / R) 2 with — Wn/v^, we obtain 



1 + ^_W2)»-1/41 



(4.3) 
(4.4) 

(4.5) 
(4.6) 

(4.7) 



[ 4 z z A ^ 

In the above we have put ^ = n + a, A = (kR) 2 / (4a), and a prime denotes the differenti- 
ation with respect to z. The general solution of fl4.7|) is given by a linear combination of 
the Whittaker functions 



W n (z) = b n M X+u /2, \u\/2(z) + C n M X+u /2, -\„\/2(z), 

where b n and arbitrary constants and M Kj ^(z) is defined by 



^ +1/2 e- 2/2 1 F 1 (/i -k+ 1/2; 2/x + 1; z) 

^+x/2 e - z/ 2 g r(2^ + i)r(/x-« + / + i/2) j 



^ r(2/x + z + i)r(/x - « + 1/2) n ' ^ 4 ' 9 ^ 

The regularity of the wave function at the origin implies that the coefficient c n of the 
singular solution M\+ u /2, -\ v \/2(z) in ( ^4.8| ) must vanish. Therefore the solution for the 
scattering region is given by 



(4.10) 



To determine the coefficients a n in ( [4.2D and b n in ( [4.10[ ), we require the continuity 
of the wave function itself as well as its derivative in the normal direction on the surface 
r = R. These conditions may be imposed on each partial wave independently to give 



a n (R) 
b n (R) 



1 inn/2 

2 



: in-w/2 



'Pn(R) , \ 

MR) J ' 

H^(kR) + e 2i5 ^ R ^H^(kR) 

M\ +u/2 , \u\/2(0i) 



(4.11) 
(4.12) 
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Here p n (R) in ( f4.11|) is denned by 

Pn (R) = kR\H^\{kR)-H { X{kR)}M x+v/2Av]/2 {a) 
-2H^\kR) {{a - 2A - v - l)M x+u/2> M/2 (a) 

+2\ + u+ l)M x+w/2+lt H/2 (a)} (4.13) 

and we have put e 2lSn ^ R ' ,a ' = p n {R)/p n {R) because it should be identified with the phase- 
shift of n-th partial wave. Thus, we obtain a solution in the scattering region 



-oo 



n=-oo 2 M \+v/2, \v\/2\0t) T 



and also in the asymptotic region 

+00 

2 



1 +00 

ifa(r,(p) = e ikrcose + - H n\ kr ) {e 2iSn{R ' a) - l} . (4.15) 



We put R — > 00 in ( 4.14| ) and ( 4.15|) by making use of the well-known asymptotic 



forms of cylindrical functions and that of the Whittaker function 

M K Jz) ~ -^r(l + 2fi)^~ 1/4 z 1/4 cos(2v^ — /itt — tt/4), (4.16) 

'IT 



for Re(/t) 3> \z\, Re(z) > in ( |4.13|) to obtain the phase-shift 



8 n (oo; a) = -(\u\ - n)n/2 . 
Then the coefficients of scattered wave in the asymptotic region is found to be 

l e inn/2 j e -<(M-n)w _ j j _ (4^7) 



When R becomes large, the coefficient b n (R) in the solution of scattering region behaves as 
b n (R) ~ e *(»-M/2> {(kR) 2 /(4a)} H/2 /V(\u\ + 1). Recalling the definition of the Whittaker 
function and a relation between the hypergeometric functions 



hm 1 F 1 (j3; T ,z/P) = oFi(r,z), (4-18) 



we observe 



^^M x+v/% H/2 (ar 2 /R 2 ) R — -^j^e'C-M^Cl + -(kr/2) 2 ) . (4.19) 
By recognizing 

( '' 7 " )H - F 1 (l + H-x 2 /4) = J wl (x) , 



rri + M 
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we finally obtain the solution for the scattering region 



*l(r> ( P)= E We*^"" 1 ^, (4-20) 

n=— oo 

as well as that for the asymptotic region 

i +00 

^ n (r ; (p) = e lkrcose + - e in(e+7T/2) {e^l"^ - l} H^(kr). (4.21) 



2 n — 



Thus we have found that in the large R limit the solution ipi(r, if) in the scattering 
region has the same form as the one that was obtained in the previous section through LS 
equation. It should be noticed, however, that in this approach the wave function ^jj(r, ip) 
describes the scattering state in the asymptotic region far from the solenoid. From ( |4.21| ) 
we find that the incident wave is given by the plane wave and that the scattered wave is 
given by the second term of (|4.21|) denoted by ^jj g(r, <p) 



1 +00 

VUs(^) = o £ e^^{e^-^-i}H^(kr). (4.22) 

n=— 00 

If we take the limit r — > 00 of ipjj g(r, <p), we again obtain the same scattering amplitude as 
in the previous section. Thus we conclude that the two approach to discuss the scattering 
problem of this system give the same result. 

Here it is better to give a comment on the relation between the methods explained here 
and in the previous section. In finding the scattering amplitude in section information 
on the scattering has been given by the asymptotic behavior of the wave function that 
corresponds to ifa(r,(p) in this section. On the other hand, in this section, ij;jj{r,ip) 
describes asymptotic behavior of the scattering state as has been shown above. The fact 
that these two methods give the same result is the consequence of the existence of the 
limit R — ► 00 in both ipr(r,(p) and iprr(r,(p) simultaneously. In other words, we may 
say that the Aharonov-Bohm scattering problem accepts the plane wave as a piece of its 
asymptotic wave function. In this regard, we are reminded of the need of a more careful 
treatment in the same analysis of the scattering by the Coulomb potential. 



5 The scattering by a solenoid with a finite radius 



As is easily seen from Q3.12 ) or (|4.20|) , the whole wave function of the scattering state 



neither vanishes nor is defined at the origin. Therefore the analyses in the preceding 
sections are unsatisfactory on this point. Fortunately the idea developed in the previous 
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section is easily generalized to the system with a solenoid of finite radius. Repeating the 
same procedure with finite r , we obtain 



^(r, tp) = £ e^~\^> \ J H (x) - ^-HS(x) \ , (5.1) 



^n(r^) = e— £ e W) i + e -.(M-^ MiiU«(4 (5.2) 

^n=-cx> [ ^H( a )J 

where a = Ax and x = kr and the wave function is assumed to vanish in the region 
< r < r . Again from the solution in the asymptotic region we can easily find the 
scattered wave 

lfe(r,¥>) = ~ E em( * W2) l + e^'-^-JgLii ^(x), (5.3) 

Z n=-oo { H \JW) 



and its asymptotic form for large x 

The first term in the scattering amplitude exactly cancels the corresponding term from 
the incident plane wave. Therefore the S'-matrix for the system is just a multiplication of 
a complex number of unit modulus: 

S n = -e-^'^ 1 ^- {\S n \ = 1), (5.5) 

on each eigenspace of the angular momentum. Hence the unitarity of the ^-matrix is 
evident. 

Let us denote a = [a] + {a} again to write 

/ r) ( OO OO ^ 

f(0) = -\j^ e ~ i[a]d | e ~ W2 E^W + eW2 E^'W) > ( 5 - 6 ) 

where A^(a) is given in terms of the Bessel and the Neumann functions by 

A+(a) = — rJ—— {cos(na/2)J n+{a} (a) - sm(na/2)N n+{a} (a)} , (5.7) 

H n+{a}( a ) 

A~(a) = — {cos(7ra/2)J n _{ a }(a) +siii(7ra/2)JV n _{ a }(a)}. (5.8) 

H n-{a}( a ) 
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The total cross section is then found to be 

4 r oo oo ^ 

*(«) = r EKWI 2 + EK»I 2 ■ (5-9) 

K ln=0 n=l J 

This result explains an interesting feature of this system: o~(a) is apparently periodic in 
a with period 2 (not 1). 

Unlike the case r = (extremely thin solenoid), the wave function is strictly subjected 
to the boundary condition ip(r Q ,(p) = even when a = integer. Thus the solenoid is 
completely impenetrable to the charged particles. Here let us consider the special case of 
a = integer. The explicit form of a is found to be 

> _ 4 J 2 (q) 8 ^ J 2 (a) 

a(eV6n) " k J » + iV » + k U + NZ(a) (5 - 10) 

for a = even integer and 

a(odd) = 4 ^ + 1 Y N ^ (5 11) 

for a = odd integer. When a tends to 0, these two formula behave in quite different 
ways. The formula ( |5.1(J| ) is nothing but a total cross section of two dimensional hard 
core scattering, thus tends to with a — > as n 2 /k{log(a/2)} 2 . This result simply means 
that in the case of a = even integer charged particles are not affected by the solenoid 
at all. On the other hand, the formula ( 5.1 1|) grows up to oo in the same limit since it 



has the Neumann function instead of the Bessel function in the numerator of each term. 
Therefore the total cross section of AB scattering for a = odd integer diverges when the 
radius of the solenoid tends to 0. This singular behavior is the common feature of the total 
cross section except the case a = even integer. If we notice that the partial cross section 
for large n approaches immediately to 4sin 2 (7ro;/2)//c even for finite a, we conclude that 
the singularity is not the consequence of putting the radius of solenoid infinitely small. 
This result implies the important fact that the vector potential can affect the charged 
particles even in the case of a = odd integer, which is the different conclusion from that 
of AB. 

Apart from the divergence of the total cross section considered above, the unitarity 



of the S'-matrix on each eigenspace of angular momentum is expected from (|5.5|) . This 
fact is also recognized from the different point of view. According to the discussion given 
in the appendix |A], the generalized optical theorem ( |A.4| ), which is rewritten in terms of 
partial wave decomposition of a scattering amplitude f(Q) = J2t=-oo e%n& fn as 

\f n \ 2 = --±=(f n + f*) , (5.12) 
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is equivalent to the unitarity of the 5-matrix. In our problem f n is given by 

fn = -^A± Ma] (a)e*™/\ (5.13) 

where the upper and the lower signs correspond ton+[a] >0orn + [ot\ < 0, respectively. 
Then the condition ( |5.12|) reads 

\At(a)\ 2 = \ {e^A±{a) + e^ 2 ^ »} , (5.14) 
and is easily verified. 

As for the system with infinitely thin solenoid, an explicit form of the S"-matrix is 
found and the operator identity S^S = SS^ = 1 can be verified directly. This is given in 
appendix ||. 



6 Results and discussions 

In this article the scattering problem, first discussed by Aharonov and Bohm [l], has 
been reconsidered. We have examined how the charged particle is scattered by the solenoid 
taking a plane wave as an incident wave in two ways ; by solving Lippmann-Schwinger 
equation ; by applying the Gordon's idea to the present situation. Furthermore the 
scattering problem by a solenoid of finite radius has been considered from the second 
viewpoint to ensure the impenetrability of the solenoid. 

We have shown that the two methods considered in this paper give the same wave 
function which represents the scattering state and that the results for the infinitely thin 
solenoid are exactly the same as was obtained by AB as far as the total wave function 
is concerned. However, it should be stressed again that the incident wave taken in this 
paper is different from that of AB to result in the disagreement in the scattering amplitude 
of this paper to that of or by the term proportional to the 5-function. It is the 
appearance of the 5-function of forward direction that guarantees the unitarity of the 
S'-matrix. On this point, our result agrees with that by Ruijsenaars [[HJ]. In appendix 
|C| we discuss the result of AB and give an explanation for the reason why we need the 
5-function in addition to the scattering amplitude given by AB. 

We also have shown that the 5-function term in the scattering amplitude makes the 
scattering cross section nonzero even in the case of a = odd integer as well as in the case 
of nonintegral a. It is better to make a comment for the case of integral a. If we regard 
a change of variables 

Tic c\ 

Mr, <p), A(x)) - (V'(r, <p), A'(x)) = (e^(r, <p), A(x) + e v ) 



16 



as a gauge transformation for integral a, it seems that there is no AB effect for that case. 
But it is not true because this causes a change in the strength of the magnetic field at the 
origin. Furthermore the new wave function is not defined at the origin, which will brake 
the single-valuedness of the wave function. We may here refer the same argument for the 
vector potential given in ||, j|. 
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A Unitarity and optical theorem in two dimensional scattering 

We provide in this appendix a note on two dimensional scattering theory for com- 
pleteness. Suppose we have two solutions, ipk(r, (p; (fo) and ipk(r, <p; (p' Q ), for a scattering 
problem corresponding to different incident beams of a same energy. They are assumed 
to have the asymptotic behavior 

Mr,<P;<Po) ~ e * rcos9 + i=e*— l*f{0) (0 = <p-<p ), (A.l) 



v T 



where the phase factor e 47r//4 has been introduced for later convenience. If we assume the 
Hamiltonian to be Hermitian, it is straightforward to obtain 

T 

d<P {ip* k (r, <p\ ip' )d r ^ k (r, ip; ip ) - ip k (r, ip; ^ Q )d r ^l{r, <p; ip' )} = (A.3) 

as a consequence of the Schrodinger equation. Taking r sufficiently large and using the 
asymptotic form of the wave functions, we immediately find 

+ J dipritP - Vb)/(V - <Po) = -/f{/(^o - Vo) + f*(<Po - ^o)} ■ (A.4) 

This is the generalized optical theorem and is nothing but the c-number version of the 
unitarity of the ^-matrix. To see this, let us define S"-matrix for the wave function given 
in ( A.l ). From the asymptotic form of the wave function 



ip k (r,(p;(p ) 




(A.5) 
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we can find a definition of operators S and / 

S=l + f, {SF){<p)=F{y) + s^£y Vo f{ V - VQ )F{ VQ ). (A.6) 
Then unitarity of the operator S reads 

Pf = -(f + P), (A.7) 

which is equivalent to (|A.4j) . 

As a special case of ( |A.4|) or ( |A.7|) , we can easily obtain the optical theorem just by 
putting (f = (p' 

a = -^2Re(/(0)). (A.8) 

B ^-matrix of the AB scattering 

Taking a plane wave 

(x\<f>(k,(po)) = J- e ifcrcose {9 = <p-<po), (B.l) 

Z7T 

as an eigenstate of the Hamiltonian(ifo) of a free particle, we obtain 

i +00 

(x\¥+\k,<p )) = — E J\u\{kry ne+i ™-^l\ (B.2) 

Z7T 



1 



n=— 00 
+00 



(X 



¥<->(*,?„)) = ^ E J\u\{kry ne ^\ (B.3) 



n=— 00 



as solutions of LS equations 

\^±\k, <p )) = \$(k, ifio)) + {E-H± ie^Vmk, <p Q )) . (B.4) 

Here we again abbreviate n + a by v. A matrix element of S operator is given in terms 
of |*(+>(fc, ¥>„)> and |*(->(fc, ip )) by 

($(p,^)|S|<%,^)> = (*(-)(p,^)|*W(g,^)) . (B.5) 

Making use of the explicit form of \^f^(k, <po)), we can easily obtain 

(tf<->(p,<p p )|tf<-%,<p,)> = E rrfr y_ <W M (pr),/ M (gr) 

n,n'=— 00 \ / " 

— in 1 Up— <p p )+in(ip— <p q )— i\v'\ir /2+in-w— i\v\ir/2 

1 1 + °° 

= S(p~q)— E e <n ^»-^ +2Wft , (B.6) 
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where 5 n = (n — \u\)n/2. It is then straightforward to see 



(^^l^l^g,^)) = j™kdkj^dy-±=5{p-k)-±=5{k-q) 



\ppk y/kq 
I +00 

x -in{tp-tp p )-2%S n in' (<p-(p q )+2iS n / 

(2tt) 2 

V / n,n'=— 00 
1 



8(p - Q)8(<Pp ~ <Pq) • ( B - 7 ) 
In the same way SS' = 1 can be verified. Therefore the ^-matrix of the AB scattering is 



unitary. By performing the sum in (|B.6|) , we can further rewrite 



(*(p,<p p )\§\*{q,<p q )) = ^S(p-q)[cosna5(9) + t S -^e^ e P- § ^ . (B.8) 



Recalling the expression (|3.18| ) for the scattering amplitude f{6) given in section we 



find a fundamental operator relation 

5 = 1+/. (B.9) 

Furthermore, if we introduce common eigenstates of Hq an d of the angular momentum 
by 

n)) = -±= 1^ d<Pe mv mk, <p)) (n = 0, ±1, ±2, . . .) , (B.10) 

V Z7T J -it 



/2tt 

the operator S is diagonalized as 



POO _ . ~ ~ 

S= kdk V \$(k,n))mk,n)\e 2iSn (B.ll) 



to convince us that the solution of LS equation assures the unitarity of S'-matrix as well 
as its commutability with H Q . 

C Result of AB and the unitarity of the ^-matrix 

The wave function of the scattering state is given by 

+00 

Mr^)= E J\v\(x)e in{e+ ^-^/\ x = kr, 6 = cp - cp , u = n + a. (C.l) 

n=— oo 

By use of the integral representation of the Bessel functions 

J v (x) = — f dte xsinht - vt [Re(x) > 0], (C.2) 



2iri Jc 



19 



we can immediately convert ( |C.1|) into its integral representation [|TTJ, [18|, 



\ p ( e ~at-ina/2 e -(l-a)t-i6+i(l+a)n /2 ~| 

^ a= 2mJc dteXSmht { i _ e-t-no+iK/2 + \ _ e -t-io+i«/2 f ' ( C - 3 ) 

for < a < 1. (See fig. |T] for the contour C.) When a has an integral part (a = [«] + {«}), 
the wave function is obtained by ip a = e -l ' Q ^ 6l+7r V{a}- We may, therefore, consider only 
the case of < a < 1. Making a change of variable, we can further rewrite ( |C.3| ) as 

= 2™ L dtC ~" FT** + 5S 1 dte "*" (C ' 4) 

where the contours C + and C_ are depicted in fig. @. On change of variable 1 1— > u = e*, 
there arises a multi-valued function w~ a in the integrand. Therefore we need to deal it 
with due care. If we recall that our solution for the scattering state, \^( + \k, ip )) in 
appendix [BL has been obtained from the LS equation, we immediately notice that we 
have only one way to deform the contour to adopt the residue theorem to the integral on 
w-plane. (See fig. |3|.) Another option for the deformation, fig. [|, obviously corresponds 
to another solution \ty(~\k, ipo)). As is seen from fig. |3|, we can make use of the residue 
theorem to the contour integration around the unit circle only when 6^0. Then we 
obtain 

ci-ri r jrr\ /*+oo p(X a )t 
^ a = e ixcose-ia(8-s g n(e)n) _ / ^ ixcosht , Q / Q x ^ 

n J-oo & — e l ° 

where 

sgn(e) = I 1 (0 < 8 < ») (C . 6) 

& v ; 1 -l (-7T < e < o) K ' 

If we interpret the modulated plane wave V'inc = e iXCOs6,_ * a ^ _sgn( - e '' )7r ' ) as an incident wave, 
the second term of ( |C.5| ) will be regarded as a scattered wave. Then we will obtain 
the Aharonov-Bohm scattering amplitude fpj${@) with the aid of the stationary phase 
approximation from 

sinvra f+°° ,. e^ 1 ^* , + 1 



7T 



jj. „«xcosht Akr—i-K /4 c / n\ 



So far the amplitude /ab(^) nas no ^ been treated in any connection with the S-matrix 
of the theory. Here it is important to note that we cannot define a S'-matrix from ( |C.5| ) 
because i/j a in ( |C.5| ) is not defined for 8 = 0. Therefore it is inappropriate to decompose 
the total wave function in the form given above for considering the relation between the 
S'-matrix and /ar(^)- ^° ^ nc ^ a definition of the S-matrix for this scattering problem, 
we need the asymptotic form of the total wave function 




e~ ix+ ™/%e + tt) + e ix ~ i7t/4 \ cos Tta6(6) + J^f AB (6) 




(C.7) 
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It should be then compared with ( |A.5|) and with discussion below in appendix [A]. For the 



present case, the S'-matrix should be defined by 

S = costal + j^B , (C.8) 
which is nothing but the result given in (|B.8[). By equating both expressions in ((J. 8) and 



that in (|B.9| ), we find 

/ = (cosvra-l)l + / AB , (C.9) 

as the relation of the two scattering amplitudes. Therefore f_£ d9\fjrQ(9)\ 2 cannot be 
interpreted as the total cross section. As a consequence, /ab(^) does no ^ obey the 
unitarity condition ( |A.7| ). Rather, it satisfies an operator relation 



/ab^ab = sin2 7ral _ cos7Ta U\B + /ab) 

because S'-matrix itself has been shown to be unitary. In terms of the amplitude itself, it 
is expressed as 

because /ab(^) satisfies 

ZabM) + /abW = o. 
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Figure 1: Schlaffi's contour of integral representaion of J v {x). C : —in + 00^ +iit + 00. 
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Figure 2: Contours C + and C_. C + : —in/2 + 00 — > +i37r/2 + 00, C_ : «7r/2 — 00 
—i3n/2 — 00 
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Figure 3: Contour on w-plane corresponding to \^( + \k, (fio)). There appears a branch cut 
on the negative Re(-u) axis. If 9 = 0, we cannot adopt the residue theorem because the 
pole is located just on the branch cut. 
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Figure 4: Another contour on w-plane corresponding to |\&( \k,(p )). The branch cut 
appears in the opposite (backward) direction for this case. 
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